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y  ABSTRACT  '  ■ 

A  simple  mathematical  model  was  developed  to  simulate 
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perfect  gases  and  the  first  law  of  thermodynamics.  Various 
input  param^t^  were  used  in  trend  analysis  to  check  the 
model.  The  output  from  the  computer  program  was  compared  to 
test  data  from  a  four  inch  bore,  open  chamber  semi-'quiescent 
diesel  engine  run  at  the  Sloan  Test  Laboratory.  This 
computer  model  was  then  modified  to  simulate  the  expansion 
stroke  of  a  newly  developed,  two  cylinder  cycle 
reciprocating  engine.  The  model  was  then  run  to  determine 
the  optimum  point  of  fuel  injection  for  the  new  engine, 
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Chanter  1 
IMTRODOCTION 


The  compression  ignition  (Cl)  engine  is  thriving  in  new 
found  popularity  amongst  automobiles,  medium-duty  and 
heavy-duty  freight  transport  trucks,  marine  propulsion  and 
auxiliary  systems,  and  various  other  industrial 
applications.  The  United  States  Department  of  Energy  (DOE) 
has  recently  predicted  that  diesel  fuel  consumption  will 
exceed  gasoline  consumption  in  this  country  by  the  year 
2000.  This  is  primarily  due  to  the  shift  to  diesel  power  in 
the  automotive  and  truck  freight  transport  industries  to 
take  advantage  of  the  high' efficiency,  high  power-to-weight 
ratio  engines.  The  present  daily  diesel  fuel  consumption  of 
the  order  of  10^  liters^^^  is  expected  to  increase  by  as 
much  as  50  percent  by  the  turn  of  the  century.  The 
increasing  fuel  consumption  rate,  coupled  with  the  constant 
concern  of  diminishing  oil  reserves,  has  prompted  renewed 
interest  in  improving  the  operating  efficiency  of  the 
conventional  compression  ignition  engine.  Any  small 
improvement  in  engine  efficiency  will  obviously  result  in  an 
enormous  savings  in  petroleum. 


The  approaches  currently  pursued  to  improve  compression 
ignition  engine  efficiency  include  increasing  the 
compression  ratio  and  the  development  of  the  "adiabatic” 
engine.  The  former  involves  turbocharging  and  the 
improvement  to  piston  ring  technology.  The  latter  approach 
concentrates  the  most  emphasis  on  insulating  the  engine. 
This  requires  the  use  of  temperature-resistant  ceramic 
cylinder  liners  for  combustion  cylinders  whose  gas  wall 
temperatures  can  be  of  the  order  or  1200  degrees  Kelvin.  In 
addition  to  these  two  approachesr  there  are  many  other 
avenues  of  research  in  progress  that  involve  improvements 
that  will  increase  the  compression  ratiOr  decrease  the  heat 
loss  from  the  engine,  or  increase  the  combustion  efficiency 
through  improved  combustion  chamber  design. 

Instead  of  improving  upon  the  conventional  compression 
ignition  engine,  a  new  cycle  engine  design  is  under 
development.  This  new  design,  proposed  by  Carmichael , 
consists  of  a  two  cylinder  cycle  which  divides  the  functions 
of  a  conventional  four-stroke  diesel  cycle  into  two  parts. 
The  new  engine  has  one  cylinder  which  compresses  the 
incoming  air  charge  and  another  cylinder  which  acts  as  the 
combustion  chamber  and  expansion  cylinder.  These  two 
cylinders  are  interconnected  by  a  regenerative  heat 
exchanger.  The  regenerator  acts  as  the  heart  of  the  new 
design.  Through  the  use  of  new  ceramic  materials,  the 
regenerator  will  act  as  a  heat  transfer  medium  by 
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tcansfering  a  portion  of  the  heat  from  the  exhaust  gases  to 
the  incoming  air  charge.  The  temperature  of  the  incoming 
air  will  be  elevated  twice.  The  first  temperature  increase 
is  due  to  the  compression  process  in  the  first  cylinder. 
This  cylinderr  in  turn,  will  transfer  its  air  charge  through 
the  ceramic  matrix  of  the  regenerative  heat  exchanger,  thus 
boosting  the  temperature  for  the  second  time.  After  passing 
through  the  regenerator,  the  incoming  air  charge  will  be  of 
sufficient  temperature  to  accomodate  spontaneous  combustion. 
With  this  high  temperature,  the  high  compression  ratio  of 
the  conventional  compression  ignition  engine  is  not  required 
to  obtain  work  from  the  cycle.  Figures  1  and  2  depict  the 
pressure  vs.  volume  and  temperature  vs.  volume  diagrams  for 
the  new  cycle  as  compared  to  a  conventional  diesel  cycle. 
The  preliminary  design  of  this  new  two  cylinder  cycle  engine 
indicates  that  an  improvement  to  thermal  efficiency  can  be 

achieved  over  the  conventional  compression  ignition  engine. 

An  important  element  of  the  engine  design  process  is 
the  capability  to  predict,  with  an  acceptable  degree  of 
accuracy,  the  energy  release  during  combustion  as  a  function 
of  time.  This  process  is  extremely  complex  in  that  it 
involves  the  injection  and  atomization  of  fuel,  the 
evaporation  and  mixing  of  the  fuel  with  the  air  charge, 
followed  by  the  various  phases  of  combustion.  The  ability 
to  accurately  predict  the  heat  release  rate  is  vital  to  the 
engine  designer  when  analyzing  a  new  engine  design. 
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This  thesis  is  an  attempt  to  assimilate  various  diesel 
engine  combustion  models  to  produce  a  simple,  yet  accurate, 
model  to  be  used  in  the  continuing  evaluation  of  the  new  two 
cylinder  cycle  reciprocating  engine.  The  proposed  model  can 
be  utilized  on  a  personal  computer  system  to  determine  the 
optimum  point  of  fuel  injection  for  the  new  engine.  The 
model  has  the  capability  to  evaluate  two  different  fuel 
types  (i  e.  iso-octane  and  propane). 
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DESCRIPTION  OF  THE  TWO  CYLINDER  CYCLE 


RECIPROCATING  ENGINE 


The  Two  Cylinder  Cycle  Reciprocating  Engine  consists  o£ 
one  compression  cylinder  and  a  pair  of  combustion/expansion 
cylinders/  see  figures  3  and  4.  The  compression  cylinder 
consists  of  an  intake  valve  and  two  exhaust  valves,  one  to 
each  expansion  cylinder.  Each  expansion  cylinder  has  its 
own  fuel  injector.  The  regenerator  cavity  contains  an 
exhaust  valve  in  addition  to  the  ceramic  matrix  regenerator. 
The  pistons  of  both  the  compression  cylinder  and  the 
expansion  cylinders  are  considered  to  be  of  simple  geometry 
with  flat  heads.  The  expansion  cylinder  piston  incorporates 
no  unique  features  to  increase  turbulence  or  swirl,  thus  it 
is  similar  to  a  direct  injection,  quiescent  chamber  diesel 
engine  cylinder.  The  five  valves  and  three  pistons  are 
actuated  by  a  camshaft  that  allows  the  compression  piston  to 
operate  at  twice  the  speed  of  an  expansion  piston.  The 
compression  cylinder  will  alternately  provide  a  compressed 
air  charge  to  each  expansion  cylinder  via  the  regenerator. 
A  typical  cycle  can  be  illustrated  by  referring  to  figures  3 
and  4. 
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step  l!  The  compression  cylinder  is  at  top  dead  center 
(TDC)  and  has  just  completed  transfering  an  air  charge  to 
cylinder  A.  Cylinder  A  is  just  starting  an  expansion  stroke 
and  cylinder  B  is  just  starting  an  exhaust  stroke.  Valve  1 
is  closed/  valve  2A  just  closed/  valve  3B  just  opened, 
valves  2B  and  3A  are  already  closed. 

Step  2;  The  compression  cylinder  is  approximately  90® 
ATDC  and  in  the  middle  of  an  air  charge  induction.  Cylinder 
A  is  still  in  the  expansion  process  and  cylinder  B  is  still 
exhausting.  Valve  1  is  open,  valves  2A,  3A,  2B  are  closed 
and  valve  3B  is  open. 

Step  3t  The  compression  cylinder  is  at  botton  dead 
center  (BDC)  and  has  completed  induction  of  an  air  charge. 
Cylinder  A  is  still  expanding  and  cylinder  B  is  exhausting 
through  the  regenerator.  Valve  1  just  closed,  valves  2A, 
3A,  2B  are  closed  and  valve  3B  is  still  open. 

Step  4 ;  The  compression  cylinder  is  approximately  90® 
BTDC  and  in  the  middle  of  compressing  the  air  charge. 
Cylinder  A  is  ending  its  expansion  stroke  and  cylinder  B  is 
completing  its  exhaust  stroke.  Valve  1  is  closed,  valves 
2A,  3A  and  2B  are  closed,  valve  3B  is  still  open. 

Step  5;  The  compression  cylinder  has  just  reached  TDC 
and  has  just  completed  an  impulse  air  charge  transfer  to 
cylinder  B  through  the  regenerator.  Cylinder  A  has  just 
commenced  its  exhaust  stroke.  Valve  1  is  closed,  valve  2B 
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has  just  closed  (it  only  opened  for  a  very  short  time  just 
before  the  compression  cylinder  reached  TDC) ,  valve  3A  just 
opened,  valves  2A  and  3B  are  closed.  (This  is  the  same  as 
step  1  except  that  cylinders  A  and  B  are  reversed.) 

Step  6 ;  The  compression  cylinder  is  approximately  90° 
ATDC  and  is  in  the  middle  of  an  air  charge  induction. 
Cylinder  A  is  still  exhausting  and  cylinder  B  is  in  the 
expansion  process.  Valve  1  is  open,  valves  2B,  3B,  2A  are 
closed  and  valve  3A  is  open.  (This  is  the  same  as  step  2 
except  that  cylinders  A  and  B  are  reversed.) 

step  7 1  The  compression  cylinder  has  just  reached  BDC 
and  has  completed  induction  of  an  air  charge.  Cylinder  A  is 
exhausting  through  the  regenerator  and  cylinder  B  is  still 
expanding.  Valve  1  just  closed,  valves  2B,  3B,  2A  are 
closed  and  valve  3A  is  open.  (This  is  the  same  as  step  3 
except  that  cylinders  A  and  B  are  reversed.) 

Step  8?  The  compression  stroke  is  approximately  90° 
BTDC  and  in  the  middle  of  compressing  an  air  charge. 
Cylinder  A  is  completing  its  exhaust  stroke  and  cylinder  B 
is  ending  its  expansion  stroke.  Valve  1  is  closed,  valve 
2B,  3B  and  2A  are  closed,  valve  3A  is  still  open.  (This  is 
the  same  as  step  4  except  that  cylinders  A  and  B  are 
reversed.) 
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step  9;  This  is  the  same  as  step  1. 


Figures  5  and  6  show  the  temperature  and  pressure  as  a 
function  of  cylinder  volume  for  a  cycle. 

The  table  below  summarizes  the  sequencing  of  the  valves 
for  a  complete  cycle  of  an  expansion  cylinder. 


Table  li  SeQuencino  of  Valves 


Valve 

1 

2h 

Ih 

Step  1: 

X 

X 

X 

X 

0 

Step  2: 

0 

X 

X 

X 

0 

Step  3: 

X 

X 

X 

X 

0 

Step  4s 

X 

X 

X 

X 

0 

Step  5: 

X 

X 

0 

X 

X  • 

Step  6: 

0 

X 

0 

X 

X 

Step  7; 

X 

X 

0 

X 

X 

Step  8: 

X 

X 

0 

X 

X 

Step  9: 

X 

X 

X 

X 

0 

where 

X  - 

Valve 

closed 

and 

0  » 

Valve 

open 

be  readily 

seen. 

the 

valve 

timing 

sequence 

rather  complex.  The  timing  sequence  must  be  such  as  to 
allow  the  impulse  transfer  of  the  air  charge  to  occur 
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without  possible  blow-down  to  the  atmosphere  or  charging  the 
wrong  cylinder.  A  shift  of  the  crank  angle  must  be 
considered  to  optimize  the  air  charge  transfer  sequence  to 
the  on-line  expansion  cylinder.  Thus,  the  valve  timing 
sequence  is  a  critical  factor  in  the  correct  and  efficient 
operation  of  this  new  engine  design  and  must  be  dealt  with 
appropriately. 
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Chapter  J 

COMBOSTION  AND  COMBOSTION  MODELING 
(An  overview) 

3.1  Description  of  Diesel  Engine  Combustion 

The  diesel  engine  combustion  process  is  exceedingly 
complex  and  not  very  well  understood.  Combustion  in  the 
diesel  engine  is  characterized  by  compression  ignition,  a 
non-uniform  fuel  and  air  distribution  in  the  combustion 
chamber,  and  a  continuous  mixing  throughout  the  period  in 
which  combustion  occurs.  Due  to  the  initial  conditions  in 
the  chamber  when  fuel  is  first  injected,  the  air  charge  in 
the  cylinder  is  of  sufficient  temperature  and  pressure  to 
support  a  chain-reaction.  However,  combustion  in  the 
compression  ignition  engine  is  governed  by  the  local 
conditions  in  each  part  of  the  charge  and  not  dependent  on 
the  spread  of  the  flame  from  one  point  to  another. 
Therefore,  the  rate  of  combustion  is  dependent  on  the  state 
and  distribution  of  the  fuel  and  upon  the  pressure  and 
temperature  within  the  cylinder. 

3.1.1  The  Phases  of  Combustion 
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Ricardo  described  the  diesel  engine  combustion  process 
as  taking  place  in  three  stages;  namely  the  delay  period,  a 
period  of  rapid  combustion,  followed  by  burning  at  a 
controlled  rate.^^^  Lyn^^^  described  the  burning  process  in 
three  slightly  different  phases.  The  first  phase  is  a 
period  of  rapid  combustion  which  lasts  for  only  three 
degrees  crank  angle.  The  second  stage  is  characterized  by  a 
decreased  rate  of  heat  release  lasting  approximately  40 
degrees  crank  angle.  The  third  period  consists  of  the  fuel 
burning  at  a  very  slow  rate  which  may  persist  through  the 
remainder  of  the  expansion  stroke. 

A  combination  of  the  descriptions  of  Ricardo  and  Lyn 
may  be  more  appropriate.  The  stages  of  combustion  could  be 
divided  into  ignition  delay,  premixed  burning,  diffusion 
controlled  combustion  and  the  tail  of  combustion. 
Figure  7  depicts  the  four  stages  of  combustion  in  a  heat 
release  diagram. 

3. 1.1.1  Ignition  Delay 

The  term  ignition  delay,  or  ignition  lag,  describes  the 
time  required  by  the  preliminary  reactions  that  occur  prior 
to  the  appearance  of  flame.  The  ignition  delay  is  broken 
down  into  a  physical  delay  and  a  chemical  delay.  The 
physical  delay  period  occurs  between  the  beginning  of  fuel 
injection  and  the  onset  of  chemical  reactions.  During  this 
period,  the  fuel  is  atomized,  vaporized,  mixed  with  air  and 
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raised  in  temperature.  This  process  is  sometimes 
collectively  referred  to  as  preparation.  The  chemical 
delay  period  immediately  follows  the  physical  delay  period 
and  terminates  at  inflammation  or  ignition.  This  period  is 
characterized  by  chemical  reactions  starting  slowly  with 
pre-flame  oxidation  of  the  fuel  followed  by  local  ignition. 

The  ignition  delay  will  vary  according  to  cylinder 
temperature,  cylinder  pressure,  the  type  of  fuel,  the 
initial  temperature  of  the  fuel,  the  characteristics  of  the 
fuel  injectors  and  the  turbulence  in  the  cylinder.  The 
physical  delay  is  small  for  light  fuels  but  can  become  the 
controlling  factor  for  heavy,  viscous  fuels.  The  physical 
delay  can  be  significantly  reduced  by  using  high  injection 
pressures  and  high  turbulence  to  expedite  the  breakup  of  the 
fuel  jet. 

Semi-empirical  relationships  have  been  developed  to 
describe  the  ignition  delay.  An  estimate  for  igniton  delay 
was  developed  by  Wolfer  in  1938:^^^ 
t  -  0.44P”^*^^exp(4650/T) 

where:  t  >  ignition  delay  in  milliseconds 

P  "  cylinder  pressure  in  atmospheres 
and  T  >  temperature  in  degrees  K  at 
ignition. 
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An  estimate  by  G.arke 


in  1970  is  quite  similar  to  that  by 
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Wolfer: 

t  ■  0.22exp(5500/T)P'^*^2”^ 

v^re:  t  •  ignitioi  delay  in  seconds 

T  •  (^linder  ten^rature  in  degrees  K 
and  P  ■  cylinder  pressure  in 


Still  another  enpirical  eiqpression  for  ignition  deley  was  devel(^>ed 
Spadaccini  and  Tevelde^^^  frcn  e3q«riments  for  NASA  in  1979  with 
diesel  fuel  in  a  stea<fy  flow  facility: 

t  -  2.43xl0“V2exp<41560/Rr) 


\d)ere  t  «  ignitioi  delay  in  seconds 
P  >  pressure  in  atmospheres 
T  «  mixture  tenf^erature  in  degrees  K 
and  R  *>  gas  ocxistant  in  atm  cmVgmole^. 

Figure  8  represents  the  effects  of  ten^rature  and  pressure  on  ignition 
delay  as  determined  from  the  estimates  by  WOlfer.  Ihe  Spadaccini  and 
Tevelde  and  Clarke  relationships  yield  somewhat  similar  results. 


When  using  ignition  delay  e:q>ressionSr  it  must  be  emphasized  that 
differences  in  engines,  fuel  properties  (especially  cetane  number) ,  fuel 
injectors  and  actual  engine  tenperatures  and  pressures  make  the 
calculaticxi  rather  eqpproximate.  lliese  formulas  are  also  very  limited  by 
their  use  of  bulk  ten^ratures,  with  no  consideration  of  local 
ccnpositicxis  or  ten^ratures.^^^^ 


3.1.1.2  gcsmixed  Burning 


In  the  premixed  burning  stage,  flame  cxxurs  at  one  or  more 
locations  and  spreads  turbulently.  Hie  rate  and  amount  of  combustion 
during  this  stage  is  directly  related  to  the  fuel  preparation  rate  and 
the  length  of  the  ignition  delay  period.  Since  this  stage  of  combustion 
is  one  of  premixed  combustion,  little  carbon  (soot)  is  produced 
resulting  in  little  radiation  heat  transfer.  However,  since  the 
combustion  rate  is  so  intense,  combustion  generated  noise  is  controlled 
by  this  stage  of  combustion.  Figure  9  depicts  premixed  burning  in  a 
cylinder. 

3.1.1.3  Diffusion  Controlled  Bumlno 

Once  the  pr^xired,  or  premixed,  fuel  has  burned,  the  combustion 
process  slows  down.  The  combustion  rate  in  this  stage  will  be  dominated 
by  the  rate  of  local  air  entrainment.  Since  the  temperature  in  the 
cylinder  is  favorable  for  i^iition  in  this  stage,  the  air/fuel  mixing 
process  will  control  the  rate  of  combustion.  This  preparation  of  the 
fuel  will  be  governed  by  the  turbulence  and  swirl  in  the  cylinder, 
lyn  estimated  that  approximately  40  percent  of  the  heat  release  from 
the  combustion  of  fuel  comes  from  this  stage.  Figure  10  shows  the 
diffusion  burning  process  in  a  cylinder. 

3.1.1.4  Conbuation  Tail  Stage 

This  last  stage  of  combustion  is  characterized  by  the  cylinder 
pressure  and  temperature  falling  as  the  expansion  process  continues. 
The  rate  of  combustion  tails  off  due  to  the  chemical  kinetic  effects  as 
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the  chemical  reaction  rate  slows.  In  this  stage,  the  reaction  rate  will 
become  the  controlling  factor  instead  of  the  air/fuel  mixing  process. 
This  stage  is  also  characterized  by  diffusion  combustion  with  a  high 
producticm  and  condxistion  of  soot  particles  with  a  resultant  high  rate 
of  radiati(»i  heat  transfer.  This  last  stage  of  combustion  can  proceed 
through  the  cooftletion  of  the  eiqsansion  stroke  and  can  contribute  i^o 
20  percent  of  the  total  heat  release.  Figure  11  represents  a  ^ical 
beat  release  rate  diagram  showing  the  four  stages  of  combustion. 

3.2  CMfeustiQO  Modeling 

The  combustion  process  is  often  considered  the  most  important 
aqpect  of  an  internal  combustion  engine,  but,  at  the  same  time,  the 
least  understood  and  most  complex.  A  mathematical  model  depicting 
combustion  would  require  good  models  of  the  fuel  system  to  include  the 
injection/fuel  pump,  the  injector  nozzles,  and  fuel  lines. 
Additioially,  models  of  fuel  atomization,  vaporization,  fuel/air  mixing, 
cylinder  air  motion,  chemical  kinetics  and  pre-mixed  and  diffusion 
mixing  would  be  required.  A  model  as  oomprdiensive  as  this  has  yet  to 
be  deveoped.  guiding  states  that  this  type  of  *combustion 

modeling  is  impossible.*  Re  justifies  this  by  pointing  out  that  the 
nunber  of  governing  restraints  and  rules  outnumber  the  degrees  of 
freedom  and,  in  addition,  the  requirements  of  low  cost,  speed  and 
accuracy  must  also  be  met.  Since  tte  conplexity  of  the  real  combustion 
process  is  so  overwhelming,  substantial  8in(>lifying  asstsoptions  must  be 
made  to  obtain  soluticxis. 
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3.2.1  Types  of  Models  and  Uses 


(13} 

Bracco  categorized  ccmbustion  models  into  three  categories 
based  on  their  uses  in  examining  different  engine  problems.  The 
categories  are  the  zero-dimensional  (or  thermocfynamic)  model,  the 
guasi-dimensional  (or  entrainment)  model,  and  the  multi-dimensional  (or 
detailed)  model. 

3. 2. 1.1  Zero-dimensional  Model^^^^ 

The  zero-dimensional  model  is  structured  around  a  thermodynamic 
analysis  of  the  engine  cylinder  contents  during  the  cycle.  The 
assuR^ions  include  one-dimensional  flew,  isentropic  adiabatic  flew 
through  nozzles  simulating  flow  past  valves,  and  unburned  mixtures  as 
mixtures  of  air,  fuel  vapor  and  residual  gases.  Specific  heats  of  the 
gas  mixture  are  modeled  using  polynomial  functions  of  tenperature. 
Con{>ression  is  assumed  to  be  adiabatic.  Ccxnbustion  assumes 
thermochemical  equilibrium  and  progressive  burning  via  mass  elements. 
The  expansion  process  assumes  thermochemical  equilibrium. 

Heat  transfer  is  modeled  using  correlations  between  the  Nusselt, 
Prandtl.  and  Reynolds  numbers  frem  heat  transfer  in  steac^  turbulent 
flew  over  flat  plates  and  pipes.  These  relationships  are  in  the  form 
of: 

NU  -  aRe^^r® 

where  a,  b,  and  c  are  obtained  frero  experimental  data  for  a  specific 
engine. 


The  combustion  process  is  generally  modeled  from  an  apparent  heat 
release  or  an  experimentally  obtained  fuel  burning  rate.  One  of  the 
most  widely  used  correlations  is  based  on  the  Wid^e  Function.  In  this 
function,  the  fuel  burned  is  expressed  as  a  fraction  of  the  total  fuel 
injected. 

PB  «  1  -  expr-KjCt) 

vrtiere  FB  »  fraction  of  fuel  burned/ total 
injected 

t  «  time  frcm  ignition 

*  shc^  factor  for  combustion  curve 
K2  •  combustion  efficiency  coefficient. 

Another  typical  function  form  is  the  cosine  function: 

x(e)  -  (1/2)  {1  -  COSTS  I  (e  -  eQ)/Ae^i} 

where  X(d)  »  mass  fraction  burned  at  crank 
angle  Q 

•  crank  angle  at  the  start  of 
combustion 

andABj^  «  bum  duration. 

There  are  numerous  other  combustion  models  that  utilize  various 
heat  release  patterns.  Some  replace  the  heat  release  curve  with  two 
straight  lines.  In  this  type  of  combustion  model,  one  line  simulates 
the  rapid  combustion  of  the  bulk  of  the  injected  fuel  and  the  other  line 
represents  the  slower  combusticxi  (tose  further  down  the  expansion 
stroke. 
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An  enopirical  model  developed  Whitehouse  and  is  based  on 

elementary  combustion  principles.  Fuel  is  assumed  to  be  prepared  for 
combustion  as  a  result  of  fuel-air  mixing.  The  reaction  rate  calculates 
the  burn  rate  in  the  premixed  stage  of  combustion.  The  preparation  rate 
becomes  governing  during  the  diffusion  burning  phase  as  the  fuel  is 
assumed  to  bum  as  rapidly  as  it  is  prepared.  (The  Whitehouse  and 
model  will  be  dealt  with  in  detail  in  a  later  chapter.) 

In  general,  thermodynamic  combustion  models  are  useful  when 
performing  a  design  trade  off  or  conparison  analysis  to  evaluate  the 
effects  of  change  in  engine  design  and  operation.  Since,  however,  the 
details  of  the  combustion  process  are  an  ii^t  to  the  model,  the  results 
can  only  indicate  what  will  transpire  if  the  engine  burns  in  the 
^lecified  manner.  These  models  cannot  address  the  feasibility  of  the 
engine  curating  in  the  prescribed  manner  because  the  details  of  the 

/TCI 

burning  process  are  not  linked  to  the  engine  design  and  operation. 
3.2.1.2  Ouasi-dimensional  Model^^^^ 

Quasi-dimmsional  models  are  also  structured  around  a  thermodynamic 
analysis  of  the  engine  cylinder  during  the  cycle.  Many  of  the  same 
assut^ims  are  utilized  to  describe  the  various  portions  of  the  process 
as  are  used  in  the  thermodynamic  model.  The  combustion  process,  on  the 
other  hand,  is  based  on  more  fundamental  physical  quantities  such  as 
turbulent  intensity,  turbulent  mixing,  jet  characteristics  in  jet  mixing 
and  the  kinetics  of  the  fuel-oxidation  process. 
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The  quasi-dimensional  models  can  be  utilized  for  the  same  purposes 
as  the  zero-dimensional  models  except  that  th^  can  new  be  used  where 
changes  in  the  combustion  process  can  be  a  dominant  factor.  The  major 
drawback  of  the  quasi-dimensional  model  is  its  inability  to  examine/  in 
detail/  the  interaction  between  fluid  flew  and  engine  geometry. 


3.2.1.3 


iSssiStr 


In  a  multi-dimensional  model/  the  governing  partial-differential 
equations  describing  conservation  of  masS/  momentum/  energy  and  specieS/ 
and  the  sub-^nodels  describing  turbulence/  chemical  kijietics/  and  etc. 
are  numerically  solved  subject  to  boundary  conditions  and  other 
restraints.  These  models  have  the  potential  for  examining  the 
interaction  between  fluid  flow  and  engine  geometry  that  is  lacking  in 
the  quasi-dimensional  model.  The  detailed  model  will  predict  engine 
performance  and  emission  characteristics  from  the  first  principles  with 
virtually  no  en^irical  relationships.  iftifortunately/  giving  the 
relevant  conservation  equations  in  three-dimeisioial/  time  dependent 
formulation/  ooqpled  with  the  state  equations  and  sub-models  leads  to  a 
oGo^ter  program  that  will  tax  ev^  the  most  capable  computer  ^stem. 
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TSB  WO  CZLSDER  CSCLE  OOBOSTICW  MODEL 


Since  the  two  cylinder  cycle  reciprocating  engine  is  a  totally  new 
concept r  combustion  modeling  can  be  even  more  difficult  than  for  a 
compression  ignition  engine.  However,  the  approach  taken  models  the 
esqpansion  cylinder  of  the  new  cycle  after  a  diesel  engine  cylinder.  Ihe 
beginning  of  the  expansion  stroke  will  simulate  a  diesel  engine  with  its 
piston  at  IDC  with  a  charge  of  air.  For  this  initial  combustion  model, 
the  air  will  be  assumed  to  be  contained  within  the  cylinder,  at 
pressure,  with  no  additional  air  added  after  ejqpansion,  as  in  the  actual 
new  engine  cycle. 

4.1  Assumptions 

I 

The  assunptions  for  this  single  zone  combustion  model  are 
essentially  those  previously  mentioned  for  the  thermodynamic  type  of 
models. 

a.  The  First  Law  of  Thermodynamics  is  used  to  establish  an  energy 
balance  to  determine  the  temperature  at  the  end  of  each  step. 

b.  The  working  fluid  is  treated  as  an  ideal  gas. 

c.  The  system  contents  are  homogeneous  and  of  uniform  temperature 
and  pressure. 
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d.  The  changes  in  gas  pretties  due  to  the  rate  of  change  of  the 
gas  confosition  are  considered  to  be  negligible. 

e.  Conbustion  is  treated  <is  a  reversible  beat  release  process. 

f.  Ccnbustion  products  are  formed  in  the  prc^rtions  according  to 
the  lav  of  perfect  combustion. 

g.  No  dissociation  of  the  products  of  combustion  occurs. 

h.  Only  four  gases  are  0(»isidered  to  be  present  and  are  varied  as 
required  for  perfect  combusti(»i. 

i.  The  incoming  air  charge  is  assumed  to  be  pure  air  plus  a 
fraction  of  the  residual  gases  remaining  in  the  cylinder. 

4.2  Thermodynamics  of  Internal  Combustion  Engines 

4.2.1  Ideal 

The  assumed  thermally  ideal  gas  obeys  the  state  equation 
EV  »  Jftr 

where  p  •  pressure 
V  ■  volume 
N  a  number  of  moles 
R  a  universal  gas  constant 
and  T  a  tea^rature. 

The  specific  gas  constant,  R,  can  be  written  in  terms  of  R  and 
the  molecular  weight  of  the  gas. 


If  the  mass  of  the  gas,  then  the  state  equation  can  be  written 

as: 

pV  ■  mPT. 


Ihe  specific  internal  energy  for  an  ideal  gas  can  be  represented  as 
a  function  of  tanperature: 


u  »  f(T) 

where  u  •  specific  internal  energy 
and  f(T)  >  functicxi  of  tenperature  dependent  on 
the  gas. 

If  the  function  f(T)  is  expressed  in  the  form  of  a  limited  power  series, 

thentl7> 

where  a^  to  a^  are  constants  which  vary 
depending  on  the  gas 

and  m  internal  energy  at  absolute  zero. 

The  specific  heat  at  constant  volume  can  be  defined  as: 

■  (dq/dT)^  -(du/dT)^. 

Thus,  following  the  same  procedures  as  for  the  internal  energy. 


The  specific  enthalpy,  h,  for  an  ideal  gas  is  given  by: 
h  ■  u  +  FT. 

It  follows  that:^^^'^®^ 

h  ■  h(T)  ■  Uq  +  y  a^T^  +  Rt. 
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At  absolute  zero, 

h  -  ho  «  Uq. 

Therefore,  for  a  perfect  gas,  the  internal  energy  varies  linearly  with 
temperature  as: 

h  ■  h^  +  Rr. 

The  ^cific  heat  at  constant  pressure,  is  defined  by: 

Cp  .  (dcj/dT)p  »  (dh/dT)p. 

For  a  perfect  gas: 

Cp  ■  +  R. 

Mow,  enthcLLpy  can  be  e^^ressed  by: 
h  -  “o  +  CpT. 

For  thermodynamic  processes  with  gases  of  const^ult  composition  and 
specific  heats  undergoing  state  changes; 

ho  -  Uo  -  0. 

Then, 

u  ■  C^T; 
h  -  CpT; 

h  -  u  -  (Cp  -  C^T  -  RT; 
and  Cp  -  Cy  -  R. 

Gas  data  are  often  given  in  terms  of  enthalpy  vice  internal  energy. 
The  convent iorud  form  is: 

h(T)/Rr  -  (h  -  hQ)/Rr 

-  ai  +  ajT  +  ajT^  +  +  a^T?*. 

and  the  internal  energy  is  expressed  as: 

u(T)/Rr  ■  (Sj^-l)  +  a2T  +  a^T^  +  a^T^  + 
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The  values  foe  the  polynomial  coefficientSr  a^  to  a^  are  provided  in 

Table  2.  Other  formulations  for  the  calculation  of  enthalpy  and 
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specific  heat  are  available  in  the  literature.  '' 


liable  2;  PPlynnmIal  fioftffifiifnte 
Bangs:  -500  r  3000  Degrees  Kelvin 


% 

00-, 

2 

3.0959 

2.73114E-03 

-7.88542E-07 

8.66002E-11 

0.0 

H2O 

3.74292 

5.65590B-04 

4.95240EH)8 

-l.ai802E-ll 

0.0 

O2 

3.25304 

6.52350E-04 

-1.49524E-07 

1.53897E-11 

0.0 

3.34435 

2.94260E-04 

1.95300E-09 

-6.57470E-12 

0.0 

^8®18 

-0.71993 

4.6426E-02 

-1.6a385E-05 

-2.67009E-09 

0.0 

1.13711 

1.45532EH)2 

-2.95876E-06 

0.0 

0.0 

4.2.2  EtQBSEtiea  Qf  as  Mixtures 

Mixtures  of  gases  obey  the  follcwing. 

a.  The  gas  mixture  as  a  whole  obeys  the  equation  of  state, 
pV  ■  HRT,  vAiere  M  is  the  total  number  of  moles  of  all  species. 

b.  The  total  pressure  of  the  mixture  is  equal  to  the  sum  of  the 
pressures  which  the  individual  components/species  exert. 

c.  The  internal  energy,  enthalpy  and  entropy  of  the  mixture  equals 
the  sum  of  the  internal  energies,  enthalpies  and  entropies  which  each 
individual  oonponent/species  would  have  if  it  separately  occupied  the 


entire  volume  of  the  mixture  at  the  same  temperature 


Thusr  for  mixtures  of  ideal  gases  the  mole  fraction  is  given  by: 


*i 


M^/M 

where  ■  moles  of  a  specie 
and  H  •  total  number  of  moles. 


Then, 

■  1.0. 

Enthalpy  is  givei  by: 

H  ■  ZM^hj^  «  MZxj^h^. 
Internal  energy  is  given  by: 

U  ■  ■  M  £ 

Specific  Heats  are  givoi  by: 

Cp  ■  Z  XiCpi 


4.2.3  The  First  Law  of  Thermodynamics 


The  emphasis  of  this  model  is  the  closed  portion  of  the  cycle. 
Thereforer  the  First  Law  of  Thermodynamic  for  closed  systons  is  simply: 
dQ  -  dH  -  dU' 

%A»re  dQ  ■  heat  energy  transfer 
df  >  work  energy  transfer 
do*  •  change  in  internal  energy. 

The  internal  energy  is  defined  by: 

O'  »  0  +  KE  PE 

tdiere  0  >  the  intrinsic  internal  energy 


-32- 


KE  >  kiiietlc  energy 
PE  ■  potential  ^lergy. 

For  a  closed  system,  we  can  assume  that  PE  «  KE  «  0. 
dQ  -  dw  «dD 

where  U  ■  M  I 

H  >  total  number  of  moles 
a  mole  fraction  of  gas  i 
Ujl^  a  specific  internal  energy  of  gas 

For  non-reacting  closed  systems,  we  can  write: 
dQ  -  dW  a  du 

v^re  dW  a  pdv  a  (  x  Xj^p)dV 
and  dU  a  Md(  £  . 

For  a  reacting  closed  ^stem,  we  can  expand  this  to: 

dQ  -  pdV  a  (Oqp  -  Ooj)  Up{T)  -  UjCT) 

where  (U^  -  -  AU_ 

op  or  o 

A  Uq  a  heat  of  reaction 
Op(T)  a  oiergy  of  products  as  a 
function  of  time 
U^(T}  a  Clergy  of  reactants  as  a 
fwction  of  temperature. 

4.3  Heat  Transfer  from  the  gas  to  the  Q^linder 


Iherefore, 


To  be  able  to  balance  the  energy  in  a  real  ^stem,  the  heat 
transfer  from  the  combustion  gas  to  the  walls  of  the  cylinder  must  be 
considered.  IWo  basic  equations  are  generally  acce^^ed  for  use  in  cycle 


calculations.  These  eu:e  the  correlations  developed  by  Annand  and 
Woschni.  Bie  relationship  by  WOschni"^®^  is  based  upon  a  forced 
convection  model. 

q/AFCjd"®  <Cj^Vp+(C2  ^^"^0^  ^P' ^ 

vAiere  C^,  €2^  ^  ^  *  constants 
A  >  area 

D  •  cylinder  bore 
p  •  pressure 

«  mean  gas  tenqperature 
■  wall  temperature 
Vp  ■  piston  velocity 
Pq  a  motoring  pressure 
p*  a  trapped  pressure 
V  a  trapped  volume 
T'  a  trapped  temperature. 

Although  WOschni's  expression  is  readily  accepted,  it  does  not 
separately  distinguish  between  convection  and  rcxUation. 

The  Annand  equation  is  also  largely  based  on  turbulent  convection, 
ttilike  the  Woschni  correlation,  Annand  claims  that  the  Reynolds  number 
is  the  major  parameter  affecting  convection.  Convection  is  the  first 
term  in  his  equation.  Ihe  second  term  in  Annand 's  equation  is  a 
radiaticxi  term  assuming  grey  bo^  radiation.  Ihus:^^^^ 
q/k  •  a(k/D)  (Re)^(Tg-i^  +  cCTg^-TW^) 

«Aiere  q  ■  heat  transfer  rate 
A  «  area 


-34- 


a,brC  s  constants 
k  ■  thermal  conductivity 
D  «  bore 

Re  a  Reynolds  Kimber  *  pVjp/ix 
P  •  density 
Vp  ■  piston  velocity 
/u  a  viscosity 

Tg  a  ten^cature  of  gas  (mean) 

OV  a  temperature  of  wall. 

The  range  of  values  for  Annand's  constants  cire: 
for  a  four  stroke  engine: 
a  a  0.26 
b  a  0.75  ±  0.15 

c  -  3.88  ±  1.39  X  10“®  J/sm^K^ 
for  a  two  stroke  engine: 
a  a  0.26 
b  a  0.64  ±  0.10 
c  a  3,03  ±  1.06  X  10~® 

Since  Annand's  equation  separates  the  convective  term  fron  the 
radiation  term»  it  is  believed  that  the  Annand  correlation  is  better 
suited  to  the  new  cycle  calculations. 

4.4  The  Combustion  Model 
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In  the  process  ot  heat  release  from  coinbustionr  both  physical  and 
effects  are  involved.  Liquid  fuel  injected  into  an  engine 
inust  be  heated,  ve^rized,  and  mixed  with  oxygen  in  the  preparation 
process  prior  to  combustion.  Once  the  fuel  is  prepared,  it  may  then 
burn  at  a  rate  controlled  by  chemical  kinetics.  It  has  been 
demonstrated  that  the  time  required  for  combustion  of  the  prepared  fuel 
is  negligible  as  centred  to  the  preparation  time. 

At  the  beginning  of  the  burning  period,  chemical  kinetics  are 
important  due  to  the  low  temperatures.  When  fuel  is  first  injected  into 
a  cylinder  of  a  diesel  engine,  the  tanperature  is  generally  such  that 
rapid  burning  will  not  occur.  Additionally,  the  heat  transferred  to  the 
incoming  fuel  causes  the  temperature  to  dre^  in  the  cylinder.  As  the 
temperature  rises  in  the  cylinder,  the  combustion  rate  rises,  thus 
increasing  the  temperature.  Ihe  heat  release  rate  continues  to  rise 
until  the  lack  of  prepared  fuel  becomes  the  controlling  factor.  When 
the  excess  prepared  fuel  is  depleted,  combustion  will  proceed  at  the 
rate  of  fuel  preparatim.  Figure  12  represents  the  effects  of 
preparation  rate  and  reaction  rate  in  premixed  burning  as  a  function  of 
crank  angle. 

4.4.1  Preparation  of  Biel 

After  injectiem,  the  fuel  is  plysically  prepared  for  combustion. 
As  mentioned  before,  this  process  involves  the  atomization,  vaporization 
and  mixing  of  the  fuel  with  air.  Ihe  rate  of  preparation  can  be  assumed 
to  be  proportional  to  the  total  surface  area  of  the  fuel  sprsy  droplets. 
If  all  the  droplets  are  assumed  to  be  of  identical  size,  then  it 
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follows: 

■  l^>rrD^ 

v^re  N^  m  Mass  of  fuel  injected 
m  Mass  of  fuel  unburned 
n  ■  number  of  fuel  drc^lets 
p  >  fuel  drc^let  density 
Oq  a  Initial  droplet  diameter 
D  a  Drc^let  dianeter. 

Die  total  area 

Area  •  •  n  n  ISH^npn)^^ 

Area  a  {6M^/npDQ^)  <6M^prt)2/3 

Area  a 

Assuming  that  the  density,  p,  atvi  initial  diameter,  D^,  are  constant, 
then  the 

Area  a 

Allowing  for  the  effect  of  oxygen  availability  on  the  mixing  of  the 
fuel,  the  preparation  rate,  FR,  can  be  written  as: 

PR  a 

where  x  a  empirical  constant 
m  a  empirical  constant 
^02  a  partial  pressure  of  oxygen 
K  a  constant. 

Die  cxmstant  K  is  a  function  of  the  characteristics  of  fuel  injection, 
air  movement  and  combustion  chamber  shape.  Typical  values  for  four 
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f. 


stroke  engine  are;^^^^ 

K  «  0.008  -  0.020 
X  -  2/3 

m  «  0.4. 

4.4.2  Reaction  q£  Biel 

Since  diesel  fuel  is  not  a  pure  substance,  it  is  in^ssible  to 
ascertain  the  exact  chemiccd  equations  involved  since  the  actued 
conpounds  in  the  fuel  are  unknown.  The  temperatures  that  are  available 
frcm  experimieits  are  only  average  cylinder  temperatures.  With  these 
approximations/estimiations,  the  equations  for  reaction  rate  are  highly 
empirical.  The  degree  of  approximation  involved  may  be  justified  due  to 
the  short  time  period  during  which  chemiical  kinetics  is  of  importance. 
Also,  the  total  fuel  that  is  burned  is  equal  to  the  amount  of  fuel  that 
is  prepared.  The  reaction  rate  equation  that  was  preposed  by  Whitehouse 
and  Way  ^7,14,21}  j^ased  on  the  Arrhenius  equation. 

R«  (K'Pq2)/<Nv^ )  y^(PR-R)dx  exp(-act/T) 

vhere  R  »  reaction  rate  per  degree  crank  angle 
K*  •  empirical  constant 
act  •  empirical  constant 
^Q2  “  partial  pressure  of  oxygen 
PR  •  preparation  rate 
N  «  engine  speed  in  rpm 
T  «  cylinder  temperature. 

The  effect  of  the  ignition  delay  period  is  incorporated  in  the  Arrhenius 
type  ejpression  e:p(-act/T).  Typical  values  of  K'  and  act  are  : 
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act  »  1,4  X  10^ 

K'  ■  1.2  X  10^®  for  two  stroke  engines 
K'  ■  65  X  10^®  for  four  stroke  engines 


The  model  was  converted  tc^  ccxnputer  code  using  TRS-80  Model  III 
Disk  Basic.  The  program  listing  is  presented  in  Appendix  B. 


In  an  effort  to  set  the  enpirical  coefficients^  the  average  value 
was  used  for  all  coefficients  that  had  a  range  of  values  for  four  stroke 
engines.  The  program  was  run  and  ccnpared  to  the  data  obtained 
Remley^^^^  in  actual  engine  testing  in  the  Sloan  Automotive  Laboratory. 
Figure  13  represents  the  pressure  versus  volume  curve  for  the  model  and 
for  the  engine  run  by  Ranley.  Appendix  A  provides  Reifications  of  the 
test  engine. 
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Chapter  5 

SEZiBCTIGN  GP  FUEL  DDBCTION  FOINT 


In  order  to  obtain  the  maximum  work  and  highest  efficiency  frcm  the 
new  two  cylinder  cycle/  the  time  of  fuel  injection  should  be  c^imized. 
To  obtain  this  optimum/  a  number  of  cycles  were  run  on  the  conputer. 

5.1  Selection  of  the  Model  Coefficients 

The  model  was  run  cissuming  the  e:^>ansion  cylinder  at  TDC  with  an 
air  charge  at  a  temperature  and  pressure  of  1090^  and  10  atmosEdieres 
while  the  engine  speed  of  850  rpn  and  air/fuel  ratio  were  held  constant. 
The  selected  fuel. was  (iso-octane)  with  a  lower  heating  value  of 

4.2  X  10^  joules/kilogram  and  a  residual  air  fraction  of  0.05. 

The  model  was  run  several  times  to  obtain  a  value  of  K  in  the 
equation: 

PR  . 

The  values  of  x  and  m  were  held  constant  at  2/3  and  0.04,  respectively/ 
as  the  values  used  for  four-stroke  diesel  engines.  When  searching  for  a 
value  of  K/  a  diffusion  combustion  period  of  70  -  120  degrees  of  crank 
angle  was  sought.  This  was  found  through  several  iterations  to  occur  at 
a  value  of  K  ■  0.012. 


The  values  of  constants  for  the  reaction  rate  equation: 

were  selected  as  the  values  for  four-stroke  diesel  engines. 

With  this  ii^t  data  and  selection  of  constantSr  the  model  yields  a 
heat  release  rate  curve  which  closely  resembles  that  described  1:^ 
Ricardo,  Lyn  and  Whitehouse  et  al,  see  figure  14.  The  premixed  burning 
phase  yields  approximately  45  percent  of  the  heat  release,  the  diffusion 
controlled  burning  phase  yields  approximately  45  percent  of  the  heat 
release  with  the  tail  of  combustion  providing  the  remaining  10  percent. 

5.2  Optimizing  Fuel  Injection 

Intuitively,  the  maximum  work  and  highest  efficiency  would  be 
expected  with  fuel  injection  and  combustion  occuring  at  TCC,  or 
ininediately  thereafter.  This,  however,  does  not  appear  to  be  the  case 
when  the  data  is  evaluated.  See  figures  15  through  21.  While  the  fuel 
injection  is  varied  from  180  (TDC)  to  205  degrees  crank  angle,  with  an 
injecticxi  period  of  20  degrees,  the  thermal  efficiency  rises.  For  fuel 
injecticxi  occurring  from  180  to  195  ciegrees,  the  tonperature  at  360 
degrees  (EDC)  is  not  sufficient  to  heat  the  regenerator  matrix  to  a 
temperature  which  will  pre-heat  the  incoming  air  cdiarge  to  1090  degrees 
Kelvin  as  specified  by  the  ii^t  <3ata.  For  fuel  injection  cx;cxiring  at 
205  degrees,  and  later,  incxxnplete  (xmbustion  will  result. 

From  this  approach,  the  optiimsa  point  of  fuel  injecition  cxrcurs  at 
200  degrees  crank  angle  for  a  fuel  injec:tion  period  of  20  degrees. 
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When  the  fuel  injection  period  is  reduced  to  10  degrees,  a  similar 
pattern  is  observed.  A  fuel  injection  point  on,  or  before,  200  degrees 
results  in  the  cylinder  gas  temperature  dropping  too  low  to  support 
sufficient  air  charge  pre-heat.  Fuel  injection  on,  or  after,  210 
degrees  results  in  incomplete  combustion.  See  figures  22  through  24. 
In  this  case,  the  optimum  point  of  fuel  injection  occurs  at  205  degrees. 
The  thermal  efficiency  for  this  case  is  higher  that  the  case  of  a  20 
degree  injection  period.  Also,  the  ^cific  fuel  consumption  is  lower 
in  the  case  of  10  degree  injection  as  compared  to  20  degree  injection. 

Through  a  similar  analysis,  the  case  of  an  air/ fuel  ratio  of  25 
yields  an  optimum  fuel  injection  point  of  195  degrees  crank  angle  for  a 
period  of  20  degrees.  For  this  air/fuel  ratio,  the  value  of  K  in  the 
preparation  rate  equation  was  selected  as  0.018  to  achieve  a  similar 
heat  release  rate  pattern. 


OGMBnS  MD  REXXMCXCKnONS 


As  can  be  readily  seen,  the  output  £ran  this  type  of  thermodynamic 
model  is  dependent  on  the  value  of  the  empirical  coefficients.  Ihe 
characteristics  of  the  heat  release  rate  curve  will  shift  as  a  function 
of  air/fuel  ratio,  tenperature,  pressure  and  engine  speed,  therefore, 
only  a  conparison  of  results  frcm  a  defined  heat  release  should  be  used 
for  qualitative  conparison  analysis. 

Since  this  ccnputer  model  was  written  for  a  personal  computer,  the 
time  required  for  one  run  is  excessively  long  for  a  detailed  conparison 
analysis.  The  run  time  for  one  run  with  five  degree  increm^ts  is 
approximately  1  hour  and  20  minutes.  A  motoring  analysis  (no 
combustion)  requires  approximately  15  minutes.  Hie  amount  of  time 
required  in  the  combusticxi  iteration  process  is  the  differoice  between 
the  two.  These  times  were  obtcdned  when  running  the  program  with  no 
remark  statements  and  elimination  of  all  unnecessary  spaces  in  the 
program.  Uidoubtedly,  the  efficiency  of  the  program  can  be  somewhat 
increased  by  utilizing  some  clever  programming  techniques.  However,  the 
use  of  a  snail  conputer  strictly  dedicated  to  a  conparison  analysis  with 
crank  angle  increments  of  one  or  two  degrees  can  occupy  the  machine  for 
an  inordinate  period  of  time. 
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6.1  Recommendations 


This  single  zone  model  allows  for  cycle  studies.  Hcwever,  a 
pcdslem  that  must  be  explored  is  the  formation  of  soot  and  gaseous 
pollutants,  niis  can  be  accomplished  by  expanding  the  model  to  a  two  or 
four  zone  model.  During  this  model  expansion,  the  effect  of 
chemical  kinetics  should  be  further  examined  to  displey  a  more  realistic 
combustion  process.  The  values  of  the  coefficients  for  the  polynomial 
expression  of  enthalpy  for  the  oti^r  products  of  ccmbustion  are  readily 
available. 

The  effects  of  heat  transfer  frcm  the  ^st^  may  be  more 
appropriately  modeled  by  the  use  of  the  widely  accepted  woschni 
correlations. The  use  of  Annand's  correlation,  however,  does  allow 
for  the  separation  of  convection  and  radiation. 

The  effects  of  mixing  of  the  air  charge  with  the  fuel  must  be 
further  explored  to  determine  the  effects  on  combustion  intensity  and 
efficiency. 

The  use  of  a  larger  computer  ^stem  would  be  most  beneficial  in  a 
comparison  analysis.  Single  runs  can  be  easily  done  on  a  personal 
computer  ^stem,  however,  maiy  runs  using  small  crank  angle  increments 
are  best,  although  more  costly,  performed  on  a  main  frame  systen  capable 
of  performing  numerous  simultaneous  calculations. 

Lastly,  to  obtain  realistic  coefficients  for  the  empirical 
constants  in  the  preparation  and  reaction  rate  equations,  experiments 
using  a  rapid  compression  machine  are  considered  appropriate.  This 


would  provide  for  realistic  data  with  minimum  cost 
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Figure  8;  Effects  of  Temperature  and  Pressure  on  Ignition  Delay 
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Figure  14:  Heat  Release  Rate  Obtained  from  Computer 
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DIESEL  ENGINE  COMBUSTION  CYCLE 


RUN 

:  CARMICHAEL  ENGINE 

INPUT 

data: 

CYLINDER 

BORE  * 

.3725  METERS 

STROKE  = 

.3725  METERS 

CONNECTING  ROD  LENGTH  =  .745  METERS 

ENGINE  SPEED  *  350  RPM 

ENGINE  COMPRESSION  RATIO  =  5 

AIR  /  FUEL  RATIO 

=  30 

TRAPPED 

PRESSURE 

=  1.01325E+06 

N/MT2 

TRAPPED 

TEMPERATURE  =  1090  DEG 

KELVIN 

RESIDUAL 

AIR  FRACTION  =  .05 

FUEL  SELECTED  FOR 

THIS  ANALYSIS 

=  CSH18  CISO 

OCTANE 

WITH  A  LOWER 

HEATING  VALUE 

=  -4.2E<-07  JQL’LES/i<( 

STOICHIOMETRIC  AIR  /  FUEL  RATIO  =  15 

.1151 

FUEL  /  AIR  EQUIVALENCE 

<PHI>  =  .503836 
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DIESEL  ENGINE  COMBUSTION  CYCLE 


run:  CARMICHAEL  ENGINE 

INPUT  data: 

CYLINDER  BORE  =  .3725  METERS 

STROKE  =  .3725  METERS 

CONNECTING  ROD  LENGTH  =  .745  METERS 

ENGINE  SPEED  =  850  RPM 

ENGINE  COMPRESSION  RATIO  =  5 

AIR  /  FUEL  RATIO  =  30 

TRAPPED  PRESSURE  =  1.01325E+06  N/MT2 

TRAPPED  TEMPERATURE  *  1090  DEG  KELVIN 

RESIDUAL  AIR  FRACTION  *  .05 

fuel:  SELECTED  FOR  THIS  ANALYSIS 


=  CSHIS  (ISO  OCTANE) 


WITH  A  LOWER  HEATING  VALUE  =  -4.2E+07  JOULES/KG 


STOICHIOMETRIC  AIR  /  FUEL  RATIO 
FUEL  /  AIR  EQUIVALENCE  (PHI)  * 
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DI 

ESEL  ENGINE  COMBUSTION  CYCLE 

RUN;  CARMICHAEL  ENGINE 

INPUT  DA 

7a: 

CYLINDER  BORE  = 

.3725  METERS 

STROKE 

*  .3725  METERS 

CONNECT 

ING  ROD  LE 

NOTH  » 

745  METERS 

ENGINE 

SPEED  *  850  RPM 

ENGINE 

COMPRESSION  RATIO  * 

5 

AIR  /  FUEL  RATIO 

=  30 

TRAPPED 

PRESSURE 

*  1.0132 

5E-06  N/Mt2 

TRAPPED 

TEMPERATURE  =  1090  DEG  KELVIN 

RESIDUAL  AIR  FRACTION  »  . 

05 

FUEL  SELECTED  FOR 

THIS  ANALYSIS  =  CSHIS  CISO 

OCTANE 

WI 

TH  A  LOWER 

HEATING 

VALUE  =  -4.2E-07  JOULES/K 

STC 

ICHIOWETRIC  AIR  /  FUEL  RATIO 

=  15.1151 

FUEL  /  AIR  EQUIVALENCE 

(PHI)  = 

.503836 
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DIESEL  ENGINE  COMBUSTION  CYCLE 
run:  CARMICHAEL  ENGINE 


INPUT  data: 

CYLINDER  BORE  =  .3725  METERS 

STROKE  =  .3725  METERS 

CONNECTING  ROD  LENGTH  =  .745  METERS 

ENGINE  SPEED  =  650  RPM 

ENGINE  COMPRESSION  RATIO  =  5 

AIR  /  FUEL  RATIO  =  30 

TRAPPED  PRESSURE  =  1.01325E'-0S  N/MT2 

TRAPPED  TEMPERATURE  =  1090  DEG  KELVIN 

RESIDUAL  AIR  FRACTION  *  .05 

FUEL  SELECTED  FOR  THIS  ANALYSIS  =  CSH16  (ISO  OCTANE) 
WITH  A  LOWER  HEATING  VALUE  =  -4.2E-07  JOULES/XG 
STOICHIOMETRIC  AIR  /  FUEL  RATIO  »  15.1151 

FUEL  /  AIR  EQUIVALENCE  (PHI)  =  .503836 
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DIESEL  ENGINE  COMBUSTION  CYCLE 


RUN;  CARMICHAEL  ENGINE 

INPUT  data; 

CYLINDER  BORE  =  .3725  METERS 

STROKE  =  .3725  METERS 

CONNECTING  ROD  LENGTH  =  .745  METERS 


ENGINE 

SPEED  =  850  RPM 

ENGINE 

COMPRESSION  RATIO  * 

5 

AIR  /  FUEL  RATIO 

*  30 

TRAPPED 

PRESSURE 

=  1.01325S*06 

N/M'r2 

TRAPPED 

TEMPERATURE  =  lOSO  DEG 

KELVIN 

RESIDUAL  AIR  FRACTION  =  . 

05 

FUEL  SELECTED  FOR 

THIS  ANALYSIS 

=  C8H18  (ISO 

OCTANE) 

WITH  A  LOWER 

KEATING 

VALUE 

*  -4.2E*07  JOULES/KG 

STOICHIOMETR 

IC  AIR  /  FUEL  RATIO 

=  15 

.  1151 

FUEL  /  AIR  EQUIVALENCE 

—A  ••.'•ra 

(PHI)  = 

.503836 
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DIESEL  ENGINE  COMBUSTION  CYCLE 


RUN;  CARMICHAEL  ENGINE 

INPUT  data; 

CYLINDER  BORE  =  .3725  METERS 

STROKE  =  .3725  METERS 

CONNECTING  ROD  LENGTH  =  .745  METERS 

ENGINE  SPEED  =  850  RPM 

ENGINE  COMPRESSION  RATIO  =  5 

AIR  /  FUEL  RATIO  =  30 

TRAPPED  PRESSURE  =  1.01325E*0S  N/Mt2 

TRAPPED  TEMPERATURE  *  1090  DEG  KELVIN 

RESIDUAL  AIR  FRACTION  =  .05 

FUEL  SELECTED  FOR  THIS  ANALYSIS  =  C8H18  (ISO  OCTANE) 
WITH  A  LOWER  HEATING  VALUE  *  -4.2E+07  JOULES/KG 
STOICHIOMETRIC  AIR  /  FUEL  RATIO  *  15.1151 

FUEL  /  AIR  EQUIVALENCE  (PHI)  =  .503836 
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DIESEL  ENGINE  COMBUSTION  CYCLE 


RUN:  CARMICHAEL  ENGINE 

INPUT  data: 

CYLINDER  BORE  = 

STROKE  = 


.3725  METERS 
.3725  METERS 


.745  METERS 


CONNECTING  ROD  LENGTH  = 

ENGINE  SPEED  =  850  RPM 

ENGINE  COMPRESSION  RATIO 
AIR  /  FUEL  RATIO  = 

TRAPPED  PRESSURE  = 

TRAPPED  TEMPERATURE 
RESIDUAL  AIR  FRACTION  «  .05 

FUEL  SELECTED  FOR  THIS  ANALYSIS 
WITH  A  LOWER  KEATING  VALUE 


30 

1.01325E-06  N/KT2 
=  1090  DEG  KELVIN 


C3H1S  (ISO  OCTANE) 
-4.2E-07  JOULES/KG 


STOICHIOMETRIC  AIR  /  FUEL  RATIO  *  15.1151 


FUEL  /  AIR  EQUIVALENCE  (PHI) 


50383G 
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-  1 
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DI 

ESEL  ENG IN 

E  COMBUSTION  CYCLE 

RUN:  CARMICHAEL  ENGINE 

INPUT  data: 

CYLINDER  BORE  = 

.3725  METERS 

STROKE 

*  .3725  K 

ETERS 

CONNECT 

ING  ROD  LENGTH  =  .745  ME 

TERS 

ENGINE 

SPEED  =  850  RPM 

ENGINE 

COMPRESSION  RATIO  =  5 

AIR  /  FUEL  RATIO 

*  30 

TRAPPED 

PRESSURE 

=  1.01325E*0G 

N/Mt2 

TRAPPED 

TEMPERATURE  =  1090  DEG 

KELVIN 

RESIDUAL  AIR  FRACTION  =  .05 

FUEL  SE 

LECTED  FOR 

THIS  ANALYSIS 

=  C8H18  (ISO 

OCTANE) 

WITH  A  LOWER 

HEATING  VALUE 

*  -4.2E-07  JOULES/KG 

STOICHIOMETRIC  AIR  /  FUEL  RATIO  =  15 

.1151 

FUEL  /  AIR  EQUIVALENCE 

(PHI)  =  .5038 
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DIESEL  ENGINE  COMBUSTION  CYCLE 


RUN;  CARMICHAEL  ENGINE 

INPUT  data: 


CYLINDER  BORE  =  .3725  METERS 

STROKE  =  .3725  METERS 

CONNECTING  ROD  LENGTH  =  .7*45  METERS 

ENGINE  SPEED  =  850  RPM 

ENGINE  COMPRESSION  RATIO  *  5 

AIR  /  FUEL  RATIO  =  30 

TRAPPED  PRESSURE  *  1.01325E*06  N/M-r2 

TRAPPED  TEMPERATURE  *  1090  DEG  KELVIN 

RESIDUAL  AIR  FRACTION  =  .05 

FUEL  SELECTED  FOR  THIS  ANALYSIS  *  C8H1S 
WITH  A  LOWER  HEATING  VALUE  =  -4.2E+ 
STOICHIOMETRIC  AIR  /  FUEL  RATIO  =  15.1151 

FUEL  /  AIR  EQUIVALENCE  (PHI)  =  .503836 


(ISO  OCTANE) 
07  JOULES/KG 
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Atgendix  A 


£)peci£ications  o£  Teat  Heaiay  Bngine 


<22) 


l^pe  of  Engine 

Bore 

Stroke 

blinder  Displaoenoit 
Connecting  Rod  Length 
GQaf>re6sion  Ratio 
Ntnber  of  Conpression  Rings 
Number  of  Oil  Rings 
Nunter  of  Inlet  Valves 
Number  of  Exhaust  Valves 
Vhlve  Diameter 
Valve  Lift 

Inlet  Valve  Timing  open/close 
Exhaust  Valve  Timing  open/close 
Diameter  of  Intake  Manifold  Pipe 
Diameter  of  Exhaust  Manifold  Pipe 


Four  Stroke 
4.0  inches 
2.5  inches 
31.41  cubic  inches 
6.25  inches 
14.3  :  1 
2 
1 
2 
2 

1.286  inches 
0.280  inches 
15°BroC/50®ABDC 
50®BBDC/15°MDC 
2.00  inches 


1.60  inches 


Data  Collected  from  Test  Run 


Engine  Speed  in  RFM 
Inlet  Pressure 
Exhaust  Pressure 
Inlet  Air  Temperature 
Air  to  Fuel  Ratio 
IMEP 

Start  of  Injection 
Ignition  Delay 
Period  of  Fuel  Injection 


1450 

13.5  inches  Hg  gage 

13.4  inches  Hg  gage 

186®  F 

25.38 

88.1  psi 

12.5°  Bmc 

5.5° 

17.5° 


^gpendix  B 

OnnfMtpr  Model 

The  conqputec  pcogrean  is  written  in  TRS-80  Model  III  Disk  Basic  and 
consists  of  a  main  program  and  nine  subroutines.  The  program  listing 
has  mmerous  remarks  statemoits  inserted  to  make  the  cdgorithm  and 
computer  code  easier  to  understand.  Since  the  program  takes  a 
considerable  length  of  time  to  run,  it  is  recomnended  that  the  remark 
statemaits  be  deleted  before  running.  Samples  of  output  are  presented 
in  figures  15  to  24. 
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These  subroutines  sre  used  to  celculste  Internsl  energy*  enthslpy  and  soles 


If  ln3ect.ion  has  occurred,  then  go  to  Combustion  Subroutine. 
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